Non-linear dispersive regime of cavity QED: The dressed dephasing model 
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Systems in the dispersive regime of cavity quantum electrodynamics (QED) are approaching the 
limits of validity of the dispersive approximation. We present a model which takes into account non- 
linear corrections to dressing of the atom by the field. We find that in presence of pure dephasing, 
photons populating the cavity act as a heat bath on the atom, inducing incoherent relaxation and 
excitation. These effects are shown to reduce the achievable signal-to-noise ratio in cavity QED 
realizations where the atom is measure indirectly through cavity transmission. 
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Cavity quantum electrodynamics (QED) studies the 
coupling of a two-level system (TLS) to one mode of 
a quantized light field. This is traditionally realized 
in systems where the TLSs arc Rydberg [l[ 0, H| or 
alkali 0, [B], Q atoms and the light field enclosed in 
a high-finesse cavity. Recently, this active field of re- 
search has attracted even more attention due to the re- 
alization of cavity QED in solid-state systems. Exam- 
ples are semiconducting 0, H| and superconducting sys- 
tems U, 12 1 where the strong coupling regime was 
achieved. Mechanical oscillators could also be used to 



reach this regime [13|, [14 1 



In this letter, we study the dispersive regime of cav- 
ity QED, where the TLS-cavity detuning is larger than 
the coupling strength. With the recent realization of 
the novel strong dispersive limit in atomic 0, 0] ; super- 
conducting [T^ and micromechanical [lrj ] systems, this 
regime offers new possibilities to study light-matter in- 
teraction at its most fundamental level. We study the 
dispersive limit by going beyond the usual linear approx- 
imation [l|, [TtJ and take into account important non- 
linear contributions. We obtain a reduced master equa- 
tion (ME) and quantum trajectory equation describing 
the TLS. We show that, due to dressing of the TLS by 
the field, pure TLS dephasing causes photons in the cav- 
ity to act as a heat bath on the TLS. Using these re- 
sults, we obtain an expression for the signal-to-noise ra- 
tio (SNR) in a homodync measurement of the TLS state. 
Contrary to initial expectations 



17, 1 



this SNR satu- 
rates with increased measurement power. This result is 
consistent with experimental results obtained in a cavity 
QED realization based on superconducting circuits (cir- 
cuit QED) [19j |. Based on this observation, we suggest 
a path to optimize the SNR of dispersive measurements. 
Our work can be applied to all cavity QED realizations 
reaching the strong coupling regime and will help in in- 
creasing the measurement fidelity [2(3| of any TLS dis- 



persively coupled to a harmonic oscillator. These results 
could, for example, be applied to superconduc ting qubit 



architectures which arc dispersively measured [211 . 

Cavity QED is described by the Jaynes-Cummings 
Hamiltonian [l[ 



H, = tiLOrCt'a 
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where and o~± are ladder operators for the photon 
field and the TLS respectively, u) r the cavity frequency, 
Lu a the TLS transition frequency and g their coupling 
strength. We are interested in the dynamics of the TLS in 
the presence of photon population of the resonator. For 
this purpose, a coherent drive on the cavity is modeled 
by the Hamiltonian H d = He m (t) (a^e~ iUmt + ae iu)mt ) , 
where e m (t) and u m denote measurement amplitude and 
frequency. 

Energy relaxation results from the coupling of the cav- 
ity a nd TLS to independent baths of harmonic oscilla- 
tors 



23] 



Hj = ih 



c\ + c 7 



dw, (2) 



where j = k or 7 represent either the cavity or the TLS 
bath with c K = a and c 7 = er_. Here, gj{u) is the cou- 
pling strength to the bath mode of frequency to. In the 
Born-Markov approximation, integrating out the baths 
leads to a TLS-cavity ME of Lindblad form with cav- 
ity decay rate k = 27r|g K (w r .)| 2 and TLS relaxation rate 
71 = T-f = 27r|g 7 (a; Q )| 2 . TLS dephasing can be mod- 
eled by coupling to a longitudinal classical fluctuating 
parameter f(t) using the Hamiltonian H v = hvf(t)o~ z , 
where v characterizes the coupling strength of the TLS 
to the fluctuations. This leads to pure dephasing of the 
TLS at a rate 7^ = 2v 2 Sf(u — * 0), with S/(w) the noise 
spectrum of f{t). 

The dispersive regime of cavity QED is realized when 
|A| = \u> a — u> r \ g. In this situation, it is use- 
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FIG. 1: (Color online) Geometrical representation of H 3 in the 
subspace £ n = {|n, g) , \n — 1, e)} with decay and dephasing 
rates for a) the bare and b) dispersive basis. The axes X D 
and Z° are the original axes which are rotated by an angle 
6 n = arctan(2Ay'n) around the Y axis such that ifP is along 
the Z direction. 



ful to move to the dispersive basis by using a unitary 
transformation to diagonalize Eq. ((T|). This is done 
using D = exp [A(7V g )(a^(T_ — acr + )] , where A(iV g ) = 
arctan (2A^/i\Q /2^fW q and N q = a T a+ |e) (e| is an oper- 
ator representing the total number of excitations. While 
this can be done exactly, here we only present the result 
to third order in the small parameter A = g/ A 



if? = DtffsD 



(3) 



« h [uj r + ( + (x + C,a)a)a z ] a) a + h(uj a + x)y, 

where 2\ = 2g 2 (l — A 2 )/A is the Stark shift per pho- 
ton. This Hamiltonian is similar to the usual dispersive 
Hamiltonian [l|, [l7| , but the third-order expansion leads 
to corrections to x, ^r, and also yields a squeezing term 
(at a) 2 of amplitude ( = —g 4 /A 3 . 

The transformation D to the dispersive basis can be 
interpreted as a rotation in the TLS-cavity space. To 
understand how this transformation will affect dissipa- 
tion, it is useful to consider £ n = {\n, g) , \n — 1, e}}, 
the n-excitation subspace. In £ n , H s takes the form 
if„ = KAa z /2 + hg-t/ndx, where a x z act as Pauli op- 
erators on £ n . As illustrated in Fig. [l] this Hamiltonian 
is diagonalized by a rotation about the Y axis of angle 
6 n = arctan(2Ay / n). The first order dispersive approxi- 
mation simply corresponds to taking 6 n fa 2\,Jn, while 
the exact transformation D corresponds to performing 
the above rotation in all subspaces £ n >o- After this trans- 
formation, the new eigenstates are entangled superposi- 
tions of TLS and cavity states. Hence, in the dispersive 
basis, the TLS is dressed by the field and acquires a pho- 
ton part. 

As illustrated in Fig. [TJ it is also important to rotate 
the system-bath Hamiltonians. We represent the TLS re- 
laxation rate 71 by an arrow pointing along the original 
axis Z D . After rotation, it acquires a component along 
the X axis, corresponding to dephasing, and the magni- 
tude of its Z component is reduced. In the same way, 



we can represent dephasing as an arrow pointing along 
X D . After rotation, it acquires a component along the Z 
axis, corresponding to upward and downward transitions 
of the TLS, and the magnitude of its X component is 
reduced. Since the rotation angle 9 n varies with n, from 
this simple picture we expect the effective TLS decay 
and dressed dephasing rates to depend on the subspace 
£ n and thus on the number of photons in the cavity. 

To obtain an effective ME for the TLS in the third- 
order dispersive approximation, we apply the transforma- 
tion D to if s , H K , Hj, H v and if^. We then integrate 
out the baths to obtain Lindblad form dissipators for the 
TLS and cavity operators (23j . Starting from this trans- 
formed ME, our goal is to obtain an effective ME for the 
TLS only. To arrive at this result, we follow Ref. [l8[ 
in tracing out the cavity degrees of freedom by first ap- 
plying a polaron-type transformation to the TLS-cavity 
ME. While in the linear dispersive approximation this re- 
sults in an exact TLS ME [18| , in the non-linear case we 
approximate the photon number operator a^a by lineariz- 
ing quantum fluctuations around the classical value. In 
this way, we obtain the following ME for the TLS density 
matrix p D in the dispersive basis 
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2 (4) 
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In this expression, tD? is the Lamb and Stark shifted 
TLS transition frequency, including contributions both 
linear [l8j and quadratic in photon number (i.e. leading 
to squeezing). The effective dephasing rate lcp a takes 
into account measurement-induced dephasing [24|, with 
corrections due to the non-linear terms. Finally, the ef- 
fective decay and excitation rates are given by 



7.1. = 7i 
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+ 7k+ lAn, Tr = 7-a«, (5) 



where j K = 2Tr\g K (ui a )\ 2 /4n cr it is the Purcell decay rate 



and 7±a 



S/ (± A)/n cr it are measurement and de- 



phasing induced excitation and relaxation rates. In these 
expressions, n cr it = A 2 /4g 2 is the critical photon num- 



ber [17| and n = P g 



photons in the cavity, with P, 
state population and n g e = 



the average number of 
the ground and excited 
r . e | 2 the photon popula- 



tion corresponding to the classical fields a g . e associated 
to the TLS ground and excited states. These satisfy 
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where A' rm = uj r — u) m + C is the rotating-frame cavity 
frequency and 5 g>e = x + 2Qn gfi the cavity pull. 

As discussed in relation to the geometrical picture of 
Fig. [TJ in the dispersive picture, the TLS relaxation rate 
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can be reduced by photon population [negative sign in 
LHS of Eq. ([5])]. Dressing of the TLS with increasing n 
renders it less sensitive to noise at u> a . However, the re- 
laxation rate is now enhanced by noise at the frequency 
A. In the same way, the TLS acquires a finite excitation 
rate jp Dressing of dephasing thus leads to the mea- 
surement photon acting as a heat bath for the TLS. We 
have carried out extensive numerical calculations which 
have shown this effective ME to accurately capture the 
dynamics of the TLS up to a photon population of the 
order of n cr it for the parameters used in Fig. [5] 

The variation of these modified rates with photon pop- 
ulation can easily be probed in cavity QED implementa- 
tions in which the TLS is measured directly dSH. In 
many realizations of cavity QED however, measurement 
of the TLS is realized indirectly by probing the signal 
transmitted through the cavity in a homodyne measure- 
ment [1, @, [lj| ■ Since the ME description does not take 
into account the result of the observation, we use quan- 
tum trajectory theory [23I HH to include this informa- 
tion and derive the evolution equation for the conditional 
state. This was done in Ref. |18| for the linear dispersive 
model and is extended here to incorporate the non-linear 
effects. For phase measurement (i.e. homodyne detection 
of the ^-quadrature (HI), we find that the state of the 
TLS conditioned on the record J(t) obeys the stochastic 
master equation (SME) 



-i^[^^(t)][^)-vT^)W t ] (7) 

where £ D - is given by Eq. ([4} and _M[c]p = (c — 
(c) t )g/2 + g(c~ (c) t )/2 is the measurement superoperator 
with (c) t = Tr[cQj (t)]. The measurement result J(t) is 
J(t) = VT^i{a z ) t + £(i), where is Gaussian white 
noise satisfying E[£(t)] = and E[£(t)£(t')} = 5{t - t'). 

The SME ([7|) is of Ito form and represents a gradual 
projective measurement of a z . The rate at which infor- 
mation about the state of the TLS comes out of the cavity 
is given by T c [(t) = r/F m cos 2 (# m ), where 
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with (3 = a e — a g , [i = a e + a g , the angles Op = arg(/3), 
6ft = arg(/u) and where ^-^6 [0,7r/2]. Furthermore, 
r] is a detection efficiency parameter. Choosing the ref- 
erence phase </> (the phase of the local oscillator) such 
that 9 m = 0, the rate of information gain about the TLS 
state is r]T m . To first order in A, this can simply be un- 
derstood by noting that r]K is the rate at which photons 
leak out of the cavity and are detected, while \/3\ 2 is the 
amount of information about the TLS state encoded in 
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FIG. 2: (Color online) a) SNR vs nii ncar /n cr it for A' rm = in 
the linear (dashed red line) and second order (full blue line) 
models. Dotted green line is Eq. ([9]). Dashed-dotted black 
line includes transients, b) Cavity pull 5 vs nu n ear/?icrit in 
the linear model (dashed red line), second order approxima- 
tion (full blue line) and exact (dotted purple line), c) SNR 
vs measurement detuning A' rm at fixed photon population 
(dashed-dotted orange line n — 10, full blue line n = 1). The 
parameters are (A, g, k, 71, y v> )/27t = (1700,170,34,0.1,0.1) 
MHz, 7±a = 2A 2 7 V and 77 = 1/80 consistant with exper- 
imental observations [l!|. These parameters correspond to 
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the photons. The second order term in r m is a correc- 
tion arising from the dressing of the TLS by the photons. 
Finally, the rate T ha (t) = nT^ sin 2 ((9 m ) in the SME is ex- 
tra non-Heiscnbcrg back-action due to the measurement. 
This rate is zero when 8 m = and is maximum when 
measuring in the opposite quadrature. 

Using the quantum trajectory equation, the TLS state 
localizes on one of the basis states as the measurement 
amplitude is increased [l8|. If this localization is faster 
than the relaxation time 7-f 1 , then the measurement re- 
sult is a faithful representation of the initial TLS state. 
If the TLS has localized in its excited state we do not ex- 
pect a single jump to the ground state with mean jump 
time 1 /71 . This is in contrast to the predictions from the 
linear dispersive model 17, 18|. Instead, telegraph noise 
due to the measurement-induced excitation rate 71 is ex- 
pected. These predictions can be experimentally tested 
by measuring the waiting time between jumps and com- 
paring to 7f and 7^. 

As is clear from the above results, increasing the mea- 
surement amplitude opens a new excitation channel 7-f 
which leads to a loss of the quantum non-demolition 
(QND) character of the dispersive measurement. This 
in turn will affect the expected SNR of a homodyne 
measurement of the field. Indeed, instead of scaling 
with power 17, [l8j], we expect the SNR to saturate 
with increasing measurement power. To demonstrate 
this in the situation where bifurcation is not important, 
we define the SNR as SNR = (SNR e + SNR ff )/2, with 
SNR e(g) = r ci /(7 Te(g ) +7ie( 9 )) for the TLS initially in 
the excited (ground) state. The rates 7| e ( g ) and 7; e (g) de- 
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pend on the TLS state through photon population which 
is evaluated from the steady-state value of Eq. (|6|) . 

Fig. [5k.) shows a plot of the SNR (solid blue line) 
as a function of measurement power, e 2 ^, and scaled to 
"linear /"-crit where ni incar = 6^/(k 2 /4 + x 2 ) is the aver- 
age photon number predicted by the linear model. For 
concreteness, the parameters chosen are taken from cir- 
cuit QED [26| , but these results apply to all cavity QED 
realizations reaching the strong coupling regime. Con- 
trary to the result obtained in the linear model (dotted 
red line), the SNR saturates with increasing power. For 
re/71 small, transients can further reduce the SNR. This 
reduction is small for the chosen parameters (dashed- 
dotted black line). The value at saturation is within a 
factor of two or three larger than observed experimen- 
tally [l^l for a transmon qubit 27 1. However, at fixed g, 
the ac-Stark shift per photon 2\ for the transmon is re- 
duced substantially by the presence of extra levels result- 
ing in a lower SNR. Considering the fact that r\ is highly 
dependent on parameters that arc difficult to extract ex- 
perimentally (e.g. number of thermal photons produced 
by the leading amplifier), the agreement with theory is 
excellent. 

The saturation of the SNR is due to the increase in 
7T + 71 ari< i t° a reduction of the cavity pull with photon 
number. Numerical simulations have shown the reduced 
dispersive ME Eq. to be very accurate up to about 
"crit ~ 25 photons. However, as shown in Fig. [2b>) cav- 
ity pull in the non-linear model (full blue line) starts 
to deviate significantly from the exact result obtained 
from diagonalization of the Jaynes-Cummings Hamilto- 
nian (dashed purple line) at about n cr ;t/2 ~ 13. Our 
predictions for the SNR therefore underestimate slightly 
the exact result in the range shown in Fig. [2k). 

As shown in Fig. [2b), the SNR is maximized by choos- 
ing a measurement frequency such that A' rm = 0. This 
value was used in Fig. [2k). For this optimal detuning and 
assuming that 0(n~| t ) terms can be neglected in Eqs. © 
and ([8]), we find 



SNR: 



Arjnn(x + 2(n) 



[re 2 /4+( X + 2Cn) 2 ][7T+7l]' 



(9) 



where we have taken <f> = 6p. This is shown as the green 
dotted line in Fig. [2k) where it is seen to be a very good 
approximation at low power. 

Using these results, it is possible to find parameters 
that maximize the SNR. For example, for a given value 
of n/n cr it one can find an optimal value of re. If 7^ + 71 
was independent of n, we find from Eq. ([9]) that the opti- 
mal value would be re opt /2 = (x + 2£n) which depends on 
power. In general, it will also depend on the relaxation 71 
and dephasing 7^ rates. In the case where 7 V ~ 71, the 
optimal re is approximately re opt , as 71 + 7^ is only weakly 
dependent on n. By contrast, when dephasing dominates 
over relaxation the optimal re increases with power. This 
is to avoid unwanted TLS mixing by limiting the photon 



population. Clearly, dephasing plays an important role 
in the reduction of the SNR and is therefore a crucial 
parameter to suppress. For superconducting qubit real- 
izations, efforts in this direction [2?| have already paid 
off [2r3 | , something which is promising in light of the cur- 
rent results. 

In summary, we introduced the dressed dephasing 
model of dispersive cavity QED. We obtained reduced 
master and quantum trajectory equations which incorpo- 
rate non- linear corrections. Dressing by the field leads 
to mixing rates for the TLS that are proportional to 
photon number and pure TLS dephasing. This dress- 
ing of dephasing reduces the QND character of disper- 
sive measurements and leads to saturation of the SNR 
with measurement power. This is in contrast to earlier 
studies and is consistant with experimental observations. 
These results apply to all physical realizations of cavity 
QED reaching the strong coupling regime and offer an 
approach to optimize the SNR (or measurement fidelity) 
in the ubiquitous TLS-harmonic oscillator system. 
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